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The first three measurements from the RHIC program were results on global observ- 
ables: charged particle multiplicity (N ch ), transverse energy (E T ) and elliptic flow (t> 2 ). 
They offer a look at the large-scale features of particle production in high-energy nuclear 
collisions, with particular insight into entropy production and collective behavior. Results 
from all of the RHIC experiments are discussed in light of data from lower energy nuclear 
collisions as well as from high-energy hadronic collisions to test our current understanding 
of the collision dynamics. 

1. Introduction 

It is hoped that Au+Au collisions at RHIC will form a large, collective state of ther- 
mally equilibrated matter. Such a state should be describable by simple quantities, such 
as energy density and pressure, which can be related to lattice or hydrodynamic calcula- 
tions incorporating a nuclear equation of state. One means to study this is to measure 
global quantities that characterize the entire event. Looking along the beam axis, we can 
attempt to understand the energy density achieved in the collision by studying the mul- 
tiplicity and Et distributions (dN c h/dr], dEx/drf). These can also be compared with pp, 
pp and pA collisions at higher and lower energies to study their scaling with the number 
of participating nucleons as well as the number of binary collisions. Azimuthal distribu- 
tions in slices of 77, studied by means of their Fourier coefficients, can be compared with 
results from lower energy nuclear collisions, where a steady increase of V2 with energy has 
been observed. This indicates that the particle emission becomes increasingly ellipsoidal, 
suggesting increased transverse pressure. Taken together, these measurements offer the 
first look at collective particle production and dynamics at RHIC. It is an open question, 
however, whether features of the initial state can survive the hadronization process. We 
review the data discussed at Quark Matter 2001 |l| to test our understanding of how the 
dynamics of the early stages of RHIC collisions can be observed in the final state. 



2. Centrality 

To understand whether particle production in nuclear collisions is fundamentally differ- 
ent than in pp or pA collisions, it is important to understand the collision geometry. The 
impact parameter of the collision determines N part , the number of nucleons that interact 
inelastically (i.e. participate), as well as N spec , the number of spectator nucleons which 
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Figure 1. Correlation between the ZDC en- 
ergy and multiplicity between 3 < \rj\ < 4.5 
(from PHOBOS). 



Figure 2. Spectrum of a single ZDC detector 
(from PHENIX). 



do not interact and continue along the beam direction. To determine N part in fixed target 
experiments, a zero-degree calorimeter (ZDC) is used to measure the forward energy and 
directly infer N part by the relationship N part + N spec = 2A. At RHIC, each experiment is 
equipped with a pair of ZDCs0] placed behind dipole magnets which sweep away charged 
particles and nuclear fragments, leaving primarily spectator neutrons. Figure [l] shows 
PHOBOS data on the relationship between charged particle production in 3 < \rj\ < 4.5 
and the neutral energy measured by the RHIC ZDCs. This data suggests that above a 
certain N part , the number of produced particles is in fact monotonically related to N part . 

Although one cannot directly measure the number of participants in a given collision, 
one can use this monotonic relationship of N part with respect to N c h or Et to relate a 
fraction of the cross section to a range in N part . This is done by means of the Glauber 
model of nuclear collisions , which allows the calculation of N part as well as the number 
of binary collisions, N co u, experienced by the colliding system as a function of the impact 
parameter. 

Of course, any attempt to measure a fraction of the total inelastic cross section is 
necessarily limited by an experiment's understanding of its trigger efficiency, since the 
systematic error increases for more peripheral events. This is complicated by the fact 
that the ZDCs are sensitive to mutual Coulomb dissociation ||, as seen in Figure |2| from 
PHENIX. However, theoretical calculations are available and current comparisons of data 
and theory |5|-[7j are in good agreement. 

Despite the difficulties in doing so, there are real benefits in estimating N part and N co u. 
Scaling pA and AA results by N part allows the comparison of results to pp data. Also, 
experimentalists can correct for the effects of fluctuations due to physics effects or detector 
acceptance which bias the measurement of the actual nuclear geometry. This facilitates 
making reliable comparisons of the RHIC experiments to each other and to SPS and AGS 
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data. 



3. Charged Particle Multiplicity 

The multiplicity of charged particles produced in heavy ion collisions arises from a 
variety of physics processes. In addition to the expected soft processes seen at lower 
energies, hard processes, nuclear shadowing, and hadronic rescattering all play a role|§. 
Each of these has an effect on the number of degrees of freedom available to the colliding 
system (i.e. the entropy). Although a value of N c h in isolation does not provide insight 
into the relative contributions of the various processes, we can attempt to disentangle 
them by systematically varying the initial conditions of the collision and comparing the 
results to pp and pA data as well as theoretical models. For this, it is useful to study the 
charged particle multiplicity per participant pair, dN c h/ dr)/\N part . 

Theoretical models of particle production in RHIC collisions broadly fall into two classes 
(see ||). The first class is based on modifying the wounded nucleon model to include hard 
processes. In these, one assumes that hard and soft processes scale with binary collisions 
and wounded nucleons, respectively: dN c h/dr) = A x N part + B x N co u. Both HIJING PJlOU 
and the eikonal model proposed by Kharzeev and Nardi (KN) [|J follow this approach. 
However, while HIJING includes additional physics effects (jet quenching and nuclear 
shadowing) which lead to a linear rise in dN c h/dr)/^N part vs. N part , the KN calculation 
only uses as input the fraction of hard processes and the PHOBOS result. This leads to 



a dependence of dN c h/dr) similar to that measured by WA98 [|13|j and WA97fT4}| which is 
well described by a simple power-law form, CN^ art . 

The other class of calculations, based on parton saturation, predict a very different 
dependence on N part . EKRT fl5||| found that a geometry-dependent saturation scale 
predicts a nearly-constant dependence of dN c h/dr)/ \N part as a function of N part . Kharzeev 
and Nardi [|J also perform a calculation based on parton saturation, including the running 
of the strong coupling constant with saturation scale, that finds dNch/dr] per participant 
to scale as \n.{Q 2 s / K 2 ) , where Q 2 is the impact-parameter dependent saturation momentum 
scale. This result agrees well with their eikonal calculation, perhaps fortuitously. 

Before comparing these predictions with experimental data, it is important to make sure 
that both theory and experiment are discussing the same quantity. To convert dN c h/dy to 
dN^/dr), it is not correct simply to apply a global correction factor, e.g. 95%. Rather, one 
must apply the proper Jacobian, dy = (3dr), to the dN c h/dy distribution for each species 
of charged particle to transform it to dN c h/dr). It is useful to expand this expression, 



dN _ dN 

dr) dy \ \m T coshy 



> (!) 



where we see that the two distributions converge at large y and uit- Thus, dN c h/dr) near 
mid-rapidity generally agrees with dN c h/dy at fixed-target experiments at CERN, but 
significant deviations (~ 15 %) arise at collider experiments, even for lighter species like 
pions. 

The first published results from RHIC were the measurement of dN c h/ dr)\ ri= Q/\N part 
as a function of energy |12[], shown in Figure |[ This provided the first opportunity to 



compare extrapolations of pp and pp collisions to AA collisions at the highest available 
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Figure 3. Energy dependence of 

dNch/drjl^^/^Npart with results from 
all four RHIC experiments at ^/snn = 130 
GeV. 
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Figure 4. PHENIX and preliminary PHO- 
BOS results for dN ch / 'dT]\ v= o/lN part 



energies. The new results show that 70% more particles are produced than at the SPS 
and 40% more than the extrapolation from pp (also shown) would predict at ^ l /sNN — 130 
GeV. This is strong evidence that particle production is not simply due to independent 
NN interactions. Instead, whatever process amplifies the production at SPS energies 
relative to pp collisions is even stronger at RHIC. At Quark Matter, the other three RHIC 
experiments presented new measurements |f?|, p!6|jr7| which are consistent with the original 
PHOBOS value within the stated systematic errors. We also show the predicted energy 
dependence from HIJING and an early EKRT calculation (both from ||10|| ). Clearly, 



the increasing theoretical uncertainties emphasize the importance of varying the collision 
energy at RHIC. 

At a fixed collision energy, one useful observable for understanding the relative role of 
hard processes is the variation of dN ch /drj with N part . The results of dN ch /drj\ v=0 / \N part 
vs. Npart is presented in Figure f| which incorporates the results from PHENIX |Tj| and 
preliminary data from PHOBOS ||, which agree very well within systematic errors. These 
measurements cut off at N part ~ 70 since the systematic error on N part grows rapidly below 
this (a situation which could be improved by colliding smaller nuclei, where even central 
events would have a smaller N part than Au+Au). Interestingly, the data disfavors both 
of the leading predictions proposed before Quark Matter. Rather than being constant or 
rising linearly with N part , as predicted by EKRT and HIJING respectively, the dependence 
on N part looks most similar to that measured at the SPS and thus is in broad agreement 
with the Kharzeev-Nardi calculations [|J. 

One can also study the full distribution of dN c h/dr], which is sensitive to all of the 
abovementioned physics effects but also gives insight into the role of longitudinal ex- 
pansion and hadronic rescattering ]Tj|. Preliminary PHOBOS data on dN c h/dT]/^N part 
scaled by \N part is shown in Figure [5| As N part increases, the distribution gets narrower. 
Indeed, forward particle production per participant actually decreases, suggesting that 
some fraction of particles are pulled back towards mid-rapidity via rescattering. Events 
generated with RQMD 2.4 |[20|| , which contains a substantial amount of such rescattering, 
qualitatively verifies this hypothesis. This seems to contradict models such as HIJING 
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Figure 5. PHOBOS dN c h/drj distributions Figure 6. Preliminary PHOBOS data on to- 
scaled by the number of participants. tal charged particle multiplicity vs. N par t 

which have no additional evolution of the system after the initial particles are produced, 
leading to no change in the particle production per participant outside of \rj\ > 3. 

By integrating the dN c h/drj distributions, PHOBOS also presented a measurement of 
the total number of charged particles as a function of N mrt flI3f , shown in Figure || Results 



from a HIJING calculation are also presented and are consistent within the 10% systematic 
errors. This is remarkable considering the variety of possible physics effects that arise out 
of the collision dynamics. 



4. Transverse Energy 

The measurement of transverse energy (Et) gives similar information as that of the 
charged multiplicity measurements, but is also sensitive to the charged particle mo- 
mentum spectrum. This gives us access to physics inaccessible to the previous mea- 
surements: the mean p?, the initial energy density (through the Bjorken formula e = 
dET/dr]/7iR 2 T pl|), and the effect of longitudinal flow (pdV work) on the evolving sys- 
tem. 

At QM2001, PHENIX presented three results |16j pertaining to the produced transverse 



energy: dE T /dN ch vs. N part , dE T /dN ch vs. y/s NN , and dE T /drj/\N part vs. N part . The 
transverse energy per charged particle is shown in Figure [F[ It appears to be constant 
over the full range of centrality, which suggests that there is no dramatic modification 
of the particle spectra occurs when moving from peripheral to central collisions. It also 
suggests that the transverse energy at mid-rapidity should scale as the charged particle 
multiplicity. This is confirmed by the measurement of total transverse energy scaled by 
N vart [fPofl . However, while the transverse energy appears to scale simply with the number 
of particles, the transverse energy per charged particle as a function of ^/snn, shown in 
Figure || shows some surprising behavior. It appears that dEx/dNch is the same at the 
AGS and at RHIC. However, the two experiments at the SPS (WA98 and NA49) suggest 
that this number might be the same, or perhaps larger, at SPS energies. This would 



seem to contradict the observations of the STAR experiment [22] at this conference that 
(Pt) increases by 20% from SPS to RHIC energies. However, other results presented at 
the conference, e.g. jet-quenching, imply that particle production at RHIC shows new 
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Figure 7. Preliminary PHENIX data and 
WA98 data on E T /N ch function of N part . 



Figure 8. Ex/N c h as a function of y/s^N- 
Along with AGS and SPS results, prelimi- 
nary PHENIX data is shown. 



features that might lead to deviations from expected behavior. 

One of the most pressing questions at this conference concerns the energy density 
achieved in RHIC collisions. Using the Bjorken formula, NA49 extracted an initial energy 
density of €bj ~ 3.2 GeV/fm 3 |23|]. Applying the same reasoning to the RHIC experiments 
would give a 70% increase. However, while the measurement of seems straightforward, 
the Bjorken formula has theoretical uncertainties which do not decrease with beam energy. 
For example, the models of particle production based on parton saturation have formation 
times which vary as the inverse of the saturation scale r ~ h/Q s ~ .1 — .2 fm. This would 
give an energy density of 16 — 20GeV/ fm 3 , as shown in M. Clearly, further theoretical 
investigations are warranted in order to reduce the ambiguity in the concept of initial 
energy density. 



5. Elliptical Flow 

One of the most striking results from the initial round of RHIC results is the magnitude 



of the elliptical flow reported by the STAR experiment p4| , which approaches the levels 
predicted by hydrodynamic calculations [p5fl . Since this initial result, the other experi- 
ments have confirmed this data and new information is now available about other aspects 



of this phenomenon |26|-pq] . 

As the two nuclei collide, the vector joining their centers defines the "reaction plane". 
While one cannot directly measure the true reaction plane, subevent analysis is used to 
determine the experimental resolution obtained by using the measured particles them- 
selves to estimate the collision angle n\l/ = arctan(X)i sin(n</>j)/ J2i cos(n0j)) [p6j] . Once 
\I/ is determined, one can extract the Fourier components v% corresponding to directed 
flow, and t; 2 corresponding to the magnitude of elliptic particle flow. For a detector with 
limited acceptance, such as the PHENIX experiment with A<p = 180°, one can use the 
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Figure 9. Energy dependence of v 2 , showing 
data ranging from Bevalac to RHIC energies. 
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Figure 10. The STAR result compared 
with PHOBOS and PHENIX. PHENIX has 
a lower cutoff of 500 MeV. 



two-particle correlation function to extract v 2 , oc 1 + J2™=i % v n cos(nA0) |28[| . In prin- 
ciple, these two methods extract the same information. However, the correlation function 
method may be sensitive to other effects, e.g. jets or HBT, that would be missed by the 
other method. 

The magnitude of v 2 varies with the energy as well as the centrality of the collision. 
The centrality dependence is controlled by the eccentricity e of the nuclear overlap region. 
It has been shown p9||30l that when the system is dilute, v 2 oc e x dN/dy, where the 
rapidity density characterizes the probability of particles to rescatter. In the limit where 
the number of scatterings becomes large, only the initial geometry is important, so v 2 oc e, 
and the proportionality constant can be predicted via hydrodynamic calculations p5]. The 
compiled data for the energy dependence of v 2 is shown in Figure |9] (from pi|). The large 
anti-flow (squeeze-out) observed in low-energy nuclear collisions is seen to change sign 
at v /J/v7v = 4GeV, turning into a continuous logarithmic rise of v 2 all the way to RHIC 
energies. It is interesting to note that v 2 at RHIC is approximately 60% higher than at 
the SPS, similar to the 70% increase in dN ch /drj already mentioned. 

Since STAR's original result |^|, three of the four RHIC experiments have measured 
the dependence of v 2 with centrality, as shown in Figure |Tl]. Both STAR and PHOBOS 
measure the event with full azimuthal acceptance and comparable event plane resolution. 
PHENIX uses the correlation function method and has a lower px cutoff of 500 MeV, 
both of which may explain why the PHENIX v 2 result is somewhat higher than the other 
two. 

The pt dependence of v 2 appears to both support the hypothesis that the central region 
of RHIC collisions shows hydrodynamic behavior, as well as suggest the appearance of jet 
quenching. STAR and preliminary PHENIX data[f2S| on v 2 (px) are shown in Figure [II] 
and are in broad agreement, at least for the basic trend. Calculations given in |32| have 
reproduced the STAR data up to pt = 2 GeV. The same calculations also predict v 2 for 
identified pions and protons and finds that protons only have non-zero flow above pt ~ 400 
MeV. Also shown in Figure |TT] are calculations of jet quenching which incorporate the 
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Figure 11. Comparison of charged parti- 
cle V2(pr) with GLV calculations including 
hydro and jet-quenching. Preliminary data 
from PHENIX and STAR are shown. 



Figure 12. vi (from PHOBOS) as a function 
of T) overlaid on top of dN^/dij for N par t = 
102 (see Figure |). 



geometry of the collision zone ||33||. They find that quenching has a dramatic effect on t> 2 
at very high-p^. Both of these are consistent with the preliminary STAR data |p6fl . 

PHOBOS presented data for t> 2 vs. t] all the way to r] = 5.4[^], shown in Figure [T^. 
Although results from three-dimensional hydrodynamic calculations are not available, 
one can still ask if v-i away from mid-rapidity scales with dN c h/dy, which is similar to 
dN c h/dr] in this region. Figure [12] shows the PHOBOS data overlaid with the most 
peripheral distribution from Figure Following a suggestion by Manly [Q, one can 
scale the multiplicity distribution to roughly match i> 2 for \rj\ > 2, to find remarkable 
agreement in the shape. The sharper peak in the flow distribution might simply be due 
to the flattening of the multiplicity distribution near r\ = suggested by the Jacobian 
shown in equation [l|. 



6. Conclusions 

We have seen that the global and flow observables measured at RHIC all show non- 
trivial collective behavior which was not predicted by any model. Indeed, the presence 
of the strong elliptic flow and the evolution of dN c h/dr] at the forward and backward 
rapidities show the effect of rescattering processes that models such as HIJING do not 
incorporate. With the upcoming 200 GeV run at RHIC, there are great opportunities 
to push our understanding even further. However, it should not be forgotten that exten- 
sive species and energy scans will be needed to make sure that we have an appropriate 
understanding of the basic features of particle production before we make detailed inter- 
pretations of the high statistics RHIC data. 
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